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Improvement, dynamical fermions, and heavy quark screening in QCD3 
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First results from simulations of improved actions for both gauge fields and staggered fermion fields in three 
dimensional QCD are presented. Fhis work provides insight into some issues of relevance to lattice theories in four 
dimensions. In particular, the renormalization of the bare lattice coupling is dramatically reduced when the tree- 
level 0(a 2 ) improved action is used. Naik improvement of the staggered fermion action produces little reduction 
in scaling violations of the rho meson mass. String breaking in the heavy quark potential in the unquenched 
theory is also clearly resolved, using Wilson loops to bound the ground state energy. 



There is a long history of lattice studies of 
three-dimensional QCD, which exhibits much of 
the physics of the four-dimensional theory, in- 
cluding linear confinement M, and spontaneous 
breakdown of a "chiral" symmetry 

QCD3 is a super-renormalizable theory and, 
with fcrmionic matter fields, only disconnected 
vacuum bubble diagrams are ultraviolet diver- 
gent. In a lattice regularization of QCD3, non- 
renormalizable operators can serve to mimick the 
effects of continuum states that are excluded by 
the cutoff [||. In addition, the bare coupling 
g 2 , (of dimension 1) and fermion mass mo must 
absorb cutoff dependent renormalizations. How- 
ever, these renormalizations vanish as a — > 0. 

The Wilson gluon action in three dimensions 
has the usual form, S Wil = -0J2 x ,n>v^ v ^ 
where in SU(2) is one-half the trace of the 
plaquette, and (3 — 4/(g 2 a). At 0(a 2 ) the tree- 
level action differs from the continuum by the op- 
erator g / t^zTW^Dj + D 2 V )F^. 

Simulations were done here using a tree-level 
0{a 2 )-accurate anisotropic gluon action H 
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^{Rpu + Rvn)] , (1) 



where i? M „ is one half the trace of the 2x1 rect- 
angle in the /1 x v plane, with £3^ = ^3 = a s /at 
(a = a s ), and £ y = a t /a s = 1,2). 

The fact that g 2 and mo are cutoff-independent 



in the continuum limit means that one obtains 
extremely simple scaling laws for physical quanti- 
ties. In particular, masses in lattice units (includ- 
ing the input bare quark masses) should satisfy 
f3(am) — constant in the continuum limit. 

Results for the string tension for the unim- 
proved action were obtained by Teper Q , and are 
shown as the solid circles in Fig. [j] (the dashed 
lines show the continuum limit estimated in ||, 
fla^/a = 1.33(1)). Results for the improved ac- 
tion obtained here (on lattices with at/a s = 1/4) 
are shown as the open squares in Fig. |l|; results 
for the static potential are shown in Fig. @. 
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Figure 1. Lattice spacing dependence of fia^/a. 

A superficially surprising aspect of the results 
is the 0(a) scaling violation evident in the Wilson 
action data for y/a/g 2 . In fact this demonstrates 
the need to renormalizc the bare coupling at finite 
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a. In four dimensions, it is known that perturba- 
tive expansions in the bare coupling are spoiled 
by large corrections; higher dimension operators 
generate 0(a a ) renormalizations of the coupling, 
from insertions in ultraviolet-divergent loop di- 
agrams. A renormalized coupling defined by a 
physical quantity absorbs the scale-independent 
renormalization, and greatly improves perturba- 
tive expansions of many quantities H. 
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Figure 2. Static potential for improved glue, with 
and without dynamical staggered fermions. 

In QCD3 higher dimension operators in the lat- 
tice action can likewise generate a large renormal- 
ization of 9q. On the other hand, loop diagrams 
in QCD3 converge as the cutoff is removed. In 
particular, one-loop effects are of 0(g 2 a), as re- 
quired by dimensional analysis. This produces 
the linear scaling violation in \fo j g\- Teper has 
done extensive calculations of 3D glueball masses 
tog |§) and showed that the ratios mc/Vcr ex- 
hibit 0(a 2 ) scaling violations, as expected, since 
\fa defines a renormalized coupling. 

The distinctive signature of the renormaliza- 
tion of the bare coupling in the 3D theory ex- 
poses an interesting feature of the lattice action. 
One might expect to see a reduction in the renor- 
malization of <7q when <Si mp is used, compared to 
simulations with the Wilson action. Remarkably, 
one finds that the renormalization is in fact al- 
most completely eliminated; the error in \fa j g$ 
at (3 = 2 is reduced from about 45% with the 
Wilson action, to less than about 5% with S[ mp . 

This is a genuinely surprising result, since all 
of the higher dimension operators present in Si mp 



should produce a leading 0(g^a) renormaliza- 
tion of the bare coupling. This includes ef- 
fects of tadpole diagrams induced by higher or- 
der terms in the link expansion (the tadpole di- 
agram in 3D is linearly divergent and, for Swili 
(l-iTriy Landau - 0.063. g 2 a + O(. g 2 a ) 2 ); the re- 
sults in Figs. |l|, || were obtained without tadpole 
renormalization. Apparently, the operator series 
in the effective action converges very rapidly in 
3D, even at scales near the lattice cutoff. 

Simulations of the static quark potential in the 
presence of dynamical Kogut-Susskind fermions 
were also done here. The staggered action in 3D 
is identical in form to the 4D action, and de- 
scribes two flavors of four-component spinors: 
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vA x )x(x) [U^(x)x(x + ft) 
Ul(x - fi)x(x - A)] + 2am V x{x)x(x). (2) 



Unquenched simulations were done on isotropic 
lattices at (3 = 2 (12 2 x 8) and at /3 = 3 (16 2 x 10). 
Hybrid molecular dynamics || were used, with 
improved glue. The bare quark mass was fixed in 
units of g 2 , with flam® = 0.3. The unquenched 
heavy-light meson mass Q was also computed at 
(3 = 3, and is shown as the dashed lines in Fig. |2|. 

String breaking is clearly demonstrated, with 
the unquenched potential approaching twice the 
heavy-light meson mass at large separations. Ex- 
cellent scaling behavior in the potential is also ob- 
served (the lattice spacing increases by 50% from 
(3 = 3 to j3 = 2), despite the fact that too is ex- 
pected to undergo a cutoff-dependent renormal- 
ization, of 0{g^a) at one-loop. Significant scaling 
violations in the unquenched potential were seen 
when the Wilson gluon action was used. 

One can roughly estimate the string break- 
ing distance i?& in physical units from this data, 
Rb/{f3a) w 2.5. If the quenched string tension in 
the 3D theory is equated with its physical value 
in 4D, y/a w 0.44 GeV, this implies Rb ~ 1.5 fm, 
which is similar to estimates in 4D 11 . The spac- 
ing at (3 



2 can also be estimated asa~ 0.3 fm. 
Only (fuzzy) Wilson loop operators were used 
in these calculations. Little evidence of string 
breaking has previously been seen using Wilson 
loops, and it has been suggested that these oper- 
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ators may have little overlap with the "broken 
string" state JhJ. If true, then working on a 
coarse lattice as was done here should prove ad- 
vantageous, as this generally improves the overlap 
of an operator with the lowest-lying state. More- 
over, this allows a more efficient probe of the cor- 
relation function at larger physical times, which 
is crucial in order to isolate the ground state; at 
(3 = 2 the unquenched effective potential V (R, T) 
reached a plateau only near T sa 4a, correspond- 
ing to a large "physical" time of about 1 fm. 

Simulations were also done with some improve- 
ment of the staggered fermion action. At 0(a 2 ) 
there are two sources of discretization error in 
the 4D action: a correction to the kinetics of 
the staggered field, and flavor-changing operators 
that must be added to cancel interactions medi- 
ated by highly virtual gluons |Lq ]. Since the 3D 
action has the same lattice symmetries JfJ, one 
expects similar tree-level 0{a 2 ) errors. 

Pion and rho meson propagators in quenched 
backgrounds were evaluated for the Kogut- 
Susskind action, and for the Naik-improved ac- 
tion, which only corrects for the kinetics of the 
fermion field (the Naik action used here is iden- 
tical in form to the 4D theory studied in |l3| ). 
Wall sources and local sinks were used, with two- 
dimensional Coulomb gauge fixing. In order to 
eliminate possible 0(g 2 a) errors due to a renor- 
malization of the bare quark mass, chiral extrapo- 
lations were done, using five quark masses (linear 
extrapolations for m 2 and m p gave good fits). 

Results for the p mass in the chiral limit are 
shown Fig. |^. Simulations were done on lattices 
with a t /a s = 1/4, and volumes from 8 2 x 16 at 
(3 = 1.5 to 24 2 x 48 at (3 = 6. The Naik action 
shows little improvement; tadpole renormaliza- 
tion also makes little difference. [Given the high 
degree of improvement of the gluon action, one 
would expect 0(a 2 ) scaling violations in m p /g 2 , 
with a sufficiently accurate chiral extrapolation. 
Unfortunately, the data is not good enough to 
discriminate between O(a) and 0(a 2 ) fits.] 

These results suggest that the dominant dis- 
cretization errors in the 3D staggered action come 
from flavor-changing interactions. This is consis- 
tent with studies of flavor symmetry breaking in 
four dimensions Il3]|. 



I am indebted to Peter Lepage for his expla- 
nation of the results reported here, from which 
much of the discussion in this paper originates. I 
also thank R. Fiebig, R. Woloshyn and M. Teper 
for fruitful conversations. 




Figure 3. Quenched rho meson mass. 
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